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Pulmonary airway closure is a potentially dangerous event that can occur in microgravity 
environments and may result in limited gas exchange for flight crew during long-term space 
flight. Repetitive airway collapse and reopening subjects the pulmonary epithelium to large, 
dynamic, and potentially injurious mechanical stresses. During ventilation at low lung volumes 
and pressures, airway instability leads to repetitive collapse and reopening. During reopening, air 
must progress through a collapsed airway, generating stresses on the airway walls, potentially 
damaging airway tissues. The normal lung can tolerate repetitive collapse and reopening. 
However, combined with insufficient or dysfunctional pulmonary surfactant, repetitive airway 
collapse and reopening produces severe lung injury. Particularly at risk is the pulmonary 
epithelium. As an important regulator of lung function and physiology, the degree of pulmonary 
epithelial damage influences the course and outcome of lung injury. In this paper we present 
experimental and computational studies to explore the hypothesis that the mechanical stresses 
associated with airway reopening inflict injury to the pulmonary epithelium. 

Experimental Investigations 

Experiments were performed in a parallel plate chamber lined with pulmonary epithelial 
cells, which was constructed as an idealized model of a collapsed segment of an airway where 
the walls are held in opposition by a viscous fluid. These experiments were conducted to 
determine whether air-liquid interfacial stresses can cause damage to epithelial cells, and to 
provide response behavior that can be correlated to the mechanical stimuli determined from 
computational investigations (below). 

In a first set of experiments, a fetal rat pulmonary epithelial cell line (CCL-149, ATCC) was 
cultured to confluence on a small (1 cm ), square region of the upper plate. The narrow channel 
was filled a model airway lining fluid. Phosphate buffered saline including 0.1 mg/mL CaCb and 
MgSC >4 (PBS) was used to model a surfactant-deficient airway lining fluid. A surfactant- 
containing airway lining fluid was approximated using Infasurf (ONY, Inc.) diluted to 1 mg/mL 
phospholipid concentration in PBS. Airway “reopening” was generated by the steady progression 
of a semi-infinite bubble of air down the length of the channel using a constant rate infusion 
pump (7 or 70 mL/min). A digital camera mounted above the channel collected sequential 
overhead images of the progressing bubble, which were used to calculate bubble velocity. Once 
removed from the apparatus, the slide was incubated with 1.2 pM Ethidium homodimer- 1 (Eth- 
1) and 1.2 pM calcein AM (Molecular Probes). For each slide, the number of injured cells was 
recorded as the average number of Eth-1 stained nuclei counted in fluorescence microscopic 
images. 

In a second experimental study, we attempted to discriminate the stress magnitude from the 
stimulus duration. To do so, the stress magnitude is modified by varying the viscosity of the 
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occlusion fluid while fixing the reopening velocity across experiments. This approach causes the 
stimulus duration to be inversely related to the magnitude of the pressure gradient. We also 
explore the mechanism for acute damage and demonstrate that repeated reopening and closure is 
shown to damage the epithelial cell layer even under conditions that would not lead to extensive 
damage from a single reopening event. 

Fluid Dynamic Simulations 

The bubble and parallel-plate flow chamber was modeled as a semi-infinite bubble 
progressing within a Hele-Shaw cell. In this model the walls were separated by a distance 2 H, 
with the semi-infinite bubble progressing in the x-direction with tip velocity U. The surface 
tension, y, was constant. The capillary number, Ca = pU/y, representing the relative importance 
of viscous to surface tension effects on the bubble determines the dynamic response of the 
system. Stokes equations, VP = pV 2 u and V • u = 0, were solved using the boundary element 
method. The interfacial stress condition applied at the air-liquid interface was I o»n \=yKn, 
where <x = -PI + p(Vu + Vu T ) was the stress tensor, h was the unit normal, and k was the 
interfacial curvature. For a given Ca , the system was simulated until a steady-state meniscus had 
developed and the stress-field and bubble geometry were determined. 

Three potentially injurious components of the stress cycle associated with bubble progression 
- the shear stress, the shear stress gradient, and the pressure gradient - were analyzed. 
Regression relationships describing the behavior of these components as a function of Ca were 
determined for very small Ca (5 x 10' 4 < Ca < 2 x 10‘ 3 ). Additionally, the thickness of the thin 
film deposited by bubble progression was estimated. Dimensionless values for the experimental 
flow conditions were extrapolated from the regression equations and redimensionalized. 

Results and Discussion 

For each condition the average number of injured cells per square centimeter was measured. 
For the saline-occluded channels, bubble progression at both velocities produced significantly 
increased numbers of injured cells when compared to the control. The slow velocity resulted in a 
66-fold increase in the number of injured cells and the fast velocity produced a 20-fold increase. 
The addition of Infasurf to the occlusion fluid reduced the number of injured cells to a level 
similar to the control. These results support the hypotheses that mechanical stresses associated 
with airway reopening injure pulmonary epithelial cells and that pulmonary surfactant in the 
normal lung protects the epithelium from injury due to airway reopening. 

The stress component that best agrees with the experimentally observed trauma is the 
maximal pressure gradient. Pressure gradients create a force imbalance on the cell membrane 
over the length of the cell. In addition, cell damage remains directly correlated with the pressure 
gradient, not the duration of stress exposure. For a low profile predominately flat cell (or region 
of a cell), the non-uniformly distributed load can depress the cell and stretch the membrane. For 
high profile cells or regions of a cell, such as the protrusion cause by the nucleus, where the 
normal forces of the cell surface are nearly opposite, the pressure gradient will pinch that region. 
The pinching can tear the membrane at the base of the protrusion or force fluid upward rupturing 
the top surface of the cell. The present study thus provides additional evidence that the 
magnitude of the pressure gradient induces cellular damage in this model of airway reopening. 
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Critical Path RoadMap 
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Gravity Effects on Ventilation 
Distribution 



Pressure 



TRANSPULMONARY PRESSURE (cmH 2 0) 



pgh 


Milic-Emili, J. Respiratory Physiology I, vol. 2, 
1974 


• (A) At FRC the lower region is less expanded, but more compliant 

so it receives larger portion of ventilation 

* (B) At RV the lower lung regions experience airway closure 


VOLUME <%) 
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Airway Closure in Microgravity 
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Airway 
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MICROGRAVITY CAUSES: 

•Regional Modification of Ventilation 
•Changes of Blood Perfusion 
•Variation in Lung Capacity 

•‘Patchy’ regions of airway collapse 

Alveoli 






Alveolar Sac 
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Related Terrestrial Syndromes 


Infant Respiratory Distress Syndrome 


>Acute Respiratory Distress Syndrome 
’Ventilator-Induced Lung Injury 
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Pulmonary Multiscale Interactions 
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Motivation 


Our goal is to determine the cause of 
reopening-induced damage, and the 
surfactant properties and airway reopening 
strategies that will allow pulmonary 
airways to be opened with minimal damage 
to the lung. 
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Stresses in Airway Reopening 
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Mechanisms of Cell 
Mechanotransduction and Damage 




3-D Surface Topography 
Influences Stress 
Distribution 



^MafriV 


Adapted from Davies, 
Physiol. Rev., 1 995 
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Cell Culture Experiments 



Lung epithelial cells were: 

• Cultured in an idealized model of small airways, 

• Exposed to a moving finger of air under reopening 
conditions, 

• Examined for cellular trauma. 
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Methods - Variable Velocity 

•Lung epithelial cells (CRL-149, ATCC) cultured to 
confluence on glass microscope slides. 

•The channel dimensions were 2.5 x 7.0 x 0.17 cm. 

•Two velocities (0.27 and 2.7 cm/s) were assessed. 


VO 


•Two occlusion fluids were assessed: 

- phosphate buffered saline (PBS) and 

- 1 mg/mL Infasurf (ONY, Inc., Buffalo, NY) in 
PBS. 


•Cellular trauma was quantified using fluorescent 
staining (Live/Dead Kit, Molecular Probes). 
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Stress Field - Rigid Channel 
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Injury by a Single Bubble Progression 

(L2 cells, Live/Dead Kit) 


Bubble Velocity 

Control 0.27 cm/s 2.7 cm/s 



Pulmonary Epithelial Cells f 

Glass Plate # (L2orA549) “r 

/ 4 - .CO 


u 

(~ 0.25 - 2.5 cm/s) 


Occlusion Fluid 

(PBS or Infasurf) 


0.16 
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Injury by a Single Bubble Progression 
(L2 cells, Live/Dead Kit) 

□ PBS ■infasurf (1 mg/mL in PBs) 



n = 5 

* p < 0.01 


Bubble Velocity [cm/s] 
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The Flow Model 


Steady Flow of a Semi-Infinite Bubble in a Channel 
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Control 0.27 2.7 

Bubble Velocity [cm/s] 


Fluid 

Speed 


T s 

(dyn/cm 2 ) 


At s 

(dyn/cm 2 ) 



AP 

(dyn/cm 2 ) 



■ 

Saline 

Slow 


15.5 


9.2 



340 



M 

Saline 

Fast 


34.3 


10.1 



170 



■ 

Infasurf 

Slow 


7.9 


3.4 



89 



■ 

Infasurf 

Fast 


17.5 


3.8 



44 



11.6 
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Predictions of Cel Normal Stresses 


G lass P late 


/Epithelial Cells 



G lass P late - 


Decreasing Velocity 


Ca = 0.05 
Ca = 0.005 
Ca = 0.0005 



(x - x men *)/H 


Film Thickness 
decreases with 
decreasing 
velocity 

The pressure 
gradient on the 
cell surface 
increases with 
decreasing 
velocity 
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Investigations of the Applied Stress 
Duration 

The variable velocity experiments 
induce stresses on cells that are not 
of constant duration. 

Hypothesis: 

The slow velocity experiments may 
induce greater damage because of 
the increased exposure time. 
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Methods - Constant Velocity 


•Human Pulmonary Epithelial Cells (A549, ATCC) cultured to 
confluence on glass microscope slides. 

•The channel dimensions were 2.5 x 7.0 x 0.17 cm. 


•A single velocity (0.34 cm/s) was applied. 

•Two viscosities were used 
^i = 8x10- 3 g/(cm s) (PBS) 

\i = 8x1 O’ 2 g/(cm s) (PBS + 14% Dextran) 

•Cellular trauma was quantified using fluorescent staining 
(Live/Dead Kit, Molecular Probes). 
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c = m£ 

Traveling-Wave Behavior r 
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Pressure Gradient, not Exposure Duration, 
Determines Damage 


U 0.34 Cm/S, Mcextran 1®M-PBS 
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Investigations of Topography 


• Our system is modeled to isolate the 
influence of epithelial topography on 
the following components of the stress 
cycle during airway reopening: 

- shear stress and shear stress gradient 

- normal stress and normal stress gradient 

J ^ bubble 

J J progression 
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Computational Model 



0 


x 

• Geometric Parameters: e=a/H A= XI H 
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Computational Model 


Boundary Element Method 



Lubrication 

Theory 

Q*/2H 

• Governing Parameter: Ca<a = — 

riju 




Normal Stress Distribution 

A7HB2. Calo.01 
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Normal Stress Distribution 

a/H 1 0.1, X/hB 2, Calo.01 





Tangential Stress Distribution 

HfflU Cal 0.01 



epithelial cell 


JA 
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Tangential Stress Distribution 

a/H = 0.1, X/H 1 2, Calo.01 
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Ca Q vs. Tangential Stress and 

Gradient s/A = 0.05 


Tangential Stress Tangential Stress Gradient 



0.001 0.01 0.1 0.001 0.01 0.1 


Ca Q = [Q*/(2H)]/(y/n) Ca Q = [Q*/(2H)]/(y/n) 
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Ca Q vs. Normal Stress Gradient 

s/aBo.05 

Normal Stress Gradient 


(dT n /dx) 



Ca Q = [Q*/(2H)]/(y/|x) 
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Equilibrium Equation of State (Infasurf) 



Infasurf Concentration [mg/mL] 


Ghadiali, S. N. and D. P. Gaver, III (2000). J. Appl. Physiol. 88: 493-506. 
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Correlation of Stress and Injury 
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Influence of Non-Equilibrium Behavior 
of Infasurf 
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Dynamic Surface Tension of Infasurf 



0.0001 0.001 0.01 0.1 1 10 


Infasurf Concentration [mg/mL] 


Ghadiali, S. N. and D. P. Gaver, III (2000). J. Appl. Physiol. 88: 493-506. 
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Conclusions 

•Combined experiments and computational 
investigations allows us to estimate the mechanical 
stresses that damage epithelial cells during reopening 

•The damaging effects from reopening are likely to be 
due to a large pressure gradient from the traveling air- 
liquid interface 

•Topological effects can increase the magnitude of 
deleterious stresses. 

Non-equilibrium surface-tension effects may increase 
damage unless concentrations are large. 
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